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The majority of the area contaminated by the Fukushima Dai-ichi nuclear power plant accident is covered by 
forest. To facilitate effective countermeasure strategies to mitigate forest contamination, we simulated the 
spatio-temporal dynamics of radiocesium deposited into Japanese forest ecosystems in 201 1 using a model 
that was developed after the Chernobyl accident in 1986. The simulation revealed that the radiocesium 
inventories in tree and soil surface organic layer components drop rapidly during the first two years after the 
fallout. Over a period of one to two years, the radiocesium is predicted to move from the tree and surface 
organic soil to the mineral soil, which eventually becomes the largest radiocesium reservoir within forest 
ecosystems. Although the uncertainty of our simulations should be considered, the results provide a basis 
for understanding and anticipating the future dynamics of radiocesium in Japanese forests following the 
Fukushima accident. 

The accident at the Fukushima Dai-ichi nuclear power plant (FDNPP) resulted from the Great East Japan 
Earthquake of 11 th March, 2011. This earthquake caused a tsunami of unprecedented magnitude that 
flooded and disabled cooling systems at the power plant, leading to catastrophic releases of radionuclides 
into the environment. The main radionuclides released into the atmosphere were iodine ( 131 I, half-life 8 days) and 
cesium ( 134 Cs, half-life 2 years, and 137 Cs, half-life 30 years), which were deposited over alarge area of Japan's land 
surface 1 " 4 . Because of their longer half-lives, tracking and forecasting the behaviour and impacts of 134 Cs and 137 Cs 
are priorities 5 ; in particular, the 137 Cs environmental contamination is expected to persist for decades. 

The majority of the contaminated area (approximately 70%) is covered by forests 6 . According to studies 
conducted in the forests contaminated by the Chernobyl accident 710 , the radiocesium deposited from the 
atmosphere onto forest canopies subsequently migrates between the structural components of the forest ecosys- 
tem, including the aboveground tree components (leaves/needles and branches), soil litter and surface organic 
layers, and the underlying mineral soil. The radiocesium in the atmosphere is first deposited on leaves/needles, 
branches, bark, and the soil surface organic layer (forest floor litter). This deposited radiocesium, primarily within 
the tree canopy, is then transferred to the soil surface via rainfall and with the shedding of dead leaves/needles and 
branches from the tree to the forest floor 11 " 15 . Over the longer term, radiocesium in the surface organic layers on 
the forest floor moves into the mineral soil, while some of the radiocesium in organic and mineral soil layers is 
absorbed by plant roots and transferred into aboveground plant organs. In broad terms, although a portion of the 
deposited radiocesium might migrate from forests via stream flow, forest ecosystems provide effective reservoirs 
for radiocesium and do not tend to discharge the radiocesium substantially in the long term. 

To develop effective countermeasures to mitigate the impacts of radioactive contamination of forests, detailed 
monitoring of the radiocesium migration and re-distribution is essential. In addition, it is useful to be able to 
predict the future dynamics of the radiocesium between forest components; modelling is the best tool for this 
prediction 16 " 19 . 

Here, we report the simulated dynamics of the radiocesium deposited in forests in Japan using a model for 
radionuclides in a forest ecosystem (RIFE1 model) 81618 . Model simulations were run using initial conditions 
based on the observed radiocesium distributions among the structural components of forests contaminated by the 
Fukushima fallout. We have simulated the radiocesium dynamics in evergreen needleleaf forests and deciduous 
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broadleaf forests over a period of 20 years from the initial contam- 
ination. Furthermore, we have combined the simulation results with 
spatial data from an airborne radionuclide survey to forecast the 
radiocesium dynamics in forest ecosystems both spatially and 
temporally. 

Results 

The time courses of 137 Cs and 134 Cs within each forest component are 
shown in Figure 1. As expected from the different half- lives of 137 Cs 
and 134 Cs, the total 137 Cs activity decreases quite slowly compared to 
that of 134 Cs (Fig. IAD for 137 Cs and Fig. 1BE for 134 Cs). The total 134 Cs 
inventory was almost equivalent to 137 Cs at the time of initial depos- 
ition (year 0) but will become less than 10% of the total initial invent- 
ory after 5 years. The total 137 Cs and 134 Cs combined inventory drops 
rapidly after the fallout and will decrease to approximately half of the 
initial fallout after approximately 10 years, primarily because of the 
radioactive decay of 134 Cs (Fig. 1CF). However, the rate at which 
the total radiocesium inventory decreases will slow after 10 years, 
when 137 Cs remains as the dominant nuclide. 

As for the total radiocesium in each component, very similar 
trends are observed in both forest types (Fig. 1CF). First, the total 
radiocesium inventory in tree components drops rapidly during the 
first two years after the fallout. The radiocesium in the soil surface 
organic layer also decreases exponentially with time. In contrast, the 
inventory within the mineral soil compartment peaks during the first 
year (Fig. IF, deciduous broadleaf forests) or second year (Fig. 1C, 
evergreen needleleaf forests) and then decreases slowly, primarily 
because of radioactive decay. Because of its shorter half-life, the 
134 Cs inventory in the mineral soil compartment starts decreasing 
approximately one year after the initial fallout (Fig. 1BE). The 137 Cs 
inventory in the mineral soil compartment peaks during the fifth/ 
sixth year (Fig. IAD). 

Coupling the modelling results with an airborne gamma activity 
survey map enabled us to predict the spatio-temporal evolution of 
the total radiocesium inventory in forest ecosystems (Fig. 2 and 
Supplementary Movie SI). The top row images in Figure 2 show 
the diminishing extent of the contaminated areas due to the physical 
decay of radiocesium over a 20 year period. For example, comparing 
the map of year 0 with the distribution after 20 years, the highly 



contaminated areas (yellow and red areas) located northwest of the 
FDNPP clearly decrease, and the blue areas, which are of relatively 
low contamination, are substantially diminished in size. The lower 
three rows show the radiocesium migration between forest compo- 
nents and the overall reduction in activity due to the physical decay of 
137 Cs and 134 Cs. The radiocesium in trees and in the soil surface 
organic layer is clearly visible at the time the fallout occurred, but 
the contamination decreases rapidly, resulting in a transfer of radio- 
activity to the mineral soil during the first two years after the fallout. 
In contrast, the mineral soil contamination becomes the largest 
radiocesium reservoir in the second year. The mineral soil radio- 
cesium inventory does not decrease significantly as a result of the 
balance between the decrease in the amount of radiocesium due to 
physical decay and its migration from the soil surface organic layer 
above the mineral soil during the early stage, and because of the 
longer half-life of 137 Cs. 

Discussion 

In European forests contaminated by the Chernobyl accident, the 
organic soil layer retained the largest portion of radiocesium and thus 
played a critical role in controlling radiocesium migration in forest 
ecosystems 717,20 . However, the data obtained from our study sites in 
the Fukushima prefecture indicate that the majority of radiocesium 
had already migrated to the mineral soil layer in 2012, thus dem- 
onstrating that the migration through the organic soil layer is more 
rapid than in European forests contaminated by fallout from 
Chernobyl. This is one of important radioecological differences 
between Fukushima and Chernobyl 7,17-20 . Our modelling study is 
not designed to reveal the underlying mechanisms of this fast migra- 
tion, and therefore further observation studies and analysis of 
physicochemical properties of forest soils in Fukushima must be 
continued. However, we hypothesize that this fast migration is a 
result of the relatively warm climate and heavy rainfall which lead 
to more rapid litter decomposition and substantially thinner organic 
soil layers than in many European forests (see below). 

Our simulations demonstrate the time courses of the total radio- 
cesium inventories and the individual radionuclides ( 137 Cs and 134 Cs) 
in evergreen needleleaf and deciduous broadleaf forests. The results 
showed that, without decontamination, the total inventory will 
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Figure 1 | Time courses of 137 Cs, 134 Cs, and total radiocesium ( 137 Cs + 134 Cs). Evergreen needleleaf forests (EN, above) and deciduous broadleaf forests 
(BD, below). The redistribution of the total initial inventory was simulated on a percentage basis. 
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Figure 2 | Predicted spatio-temporal distributions of the radiocesium deposited into forests. The results at 0, 1, 3, 5, 10, and 20 years after the 
accident are shown from left to right. The upper row shows the total radiocesium in forest ecosystems; the second, third, and fourth rows show the 
radiocesium in tree (Tree), soil surface organic layer (Org.L), and mineral soil (Soil) components, respectively. The black triangle indicates the location of 
Fukushima Dai-ichi nuclear power plant (FDNPP) . White areas are not covered by forests. The spatial radiocesium distribution map we obtained did not 
include the area near the FDNPP. The maps were created using the Generic Mapping Tools (GMT) 29 . 



decrease to approximately 50% of the total initial inventory in 5-10 
years after the fallout due to the radioactive decay of 134 Cs and that, 
after 5-10 years, the mineral soil compartment will retain most of the 
radiocesium (mostly 137 Cs at > 5 years). Furthermore, the simulation 
results provide guidance for the design of effective forest ecosystem 
decontamination strategies 21 . First, despite the different radiocesium 
distributions in the two forest types at the time of fallout, the overall 
trends of the contamination inventories were the same in both ever- 
green needleleaf forests and deciduous broadleaf forests. Second, it is 
clear that removing tree components is not an effective means of 
removing the radiocesium from the ecosystems unless it is performed 
shortly after the fallout. Our simulation results suggest that the radio- 
cesium within the tree component decreases rapidly within the first 
two years. Third, removing the soil surface organic layer would be 
effective during a period of approximately two years after the fallout, 
but the effectiveness of this countermeasure decreases rapidly after 
this time. 

Another conclusion that can be drawn from our simulation is that 
the risk of radiocesium migration from forest ecosystems to adjacent 
land areas will decrease rapidly during the first two years because 
migration occurs from the soil surface organic layer into the mineral 
soil compartment, which is much less likely to be subject to lateral 
transport by wind or water. The soil surface organic layer is more 
susceptible to physical redistribution by erosion than the mineral 
soil, and furthermore, the radiocesium in mineral soils can be 
strongly fixed in mineral particles and is, thus, chemically less 
mobile 18 ' 19 . 



The forests in Japan experience high mean annual precipitation 
(approximately 1500 mm; area of -— 137, 134 Cs 50 kBq m~ 2 in 
September 2011) with occasional heavy rainfall from typhoons and 
the monsoon. This precipitation, coupled with a warm climate (mean 
annual temperature, approximately 10°C), will most likely increase 
the rate of radiocesium circulation within forest ecosystems. In par- 
ticular, it is possible that high and intense rainfall and rapid decom- 
position of organic matter within the soil surface organic layer will 
combine to accelerate radiocesium migration from the soil surface 
organic layer to the mineral soil. 

A possibly accelerated rate of radiocesium migration in Japanese 
soil surface organic layers was included in our simulations by assum- 
ing reduced thicknesses of the layers compared with those reported 
in the Chernobyl studies, and also by shorter residence half-lives for 
radiocesium within the organic soil layer (k 6 ) (see Methods and 
Supplementary Table SI). These results imply that forest decontam- 
ination by removing contaminated aboveground components (e.g., 
leaves, tree boles, branches, and soil surface organic layers) should be 
performed more rapidly than expected based on the experience of 
Chernobyl 21 . In the worst case, it might already be too late to remove 
the bulk of the forest contamination by removing the aboveground 
components, and probably the mineral soil must be removed to 
achieve effective decontamination. However, removing the mineral 
soil will incur significantly greater cost than removing the above- 
ground components. Furthermore, it should be noted that large scale 
removal of trees and soils may increase soil erosion and landslides, 
and will also result in reduction of biodiversity and degradation of 
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soil nutrient and water retention functions. Therefore, the feasibility 
of removal of contaminated materials from large areas needs to be 
evaluated by cost-benefit analysis. More optimistically, the rapid 
radiocesium migration to mineral soil components may indicate 
an equal reduction in the physico-chemical mobility of radiocesium 
in forest ecosystems. The radiocesium in mineral soil components, 
especially clay minerals, is the most stable and recalcitrant form of 
radiocesium in the environment. 

In this study, the time courses of radiocesium in forests were 
simulated using the RIFE1 model. Although the uncertainties in 
our simulations should be considered, our results support the idea 
that selective countermeasure strategies can be developed to reduce 
forest contamination. There is a wide range of different counter- 
measure options which have previously been proposed 7 ' 22 " 25 and each 
of these should be examined and tested carefully before application. 
Specifically, the potential for reduction of radiation dose, the feas- 
ibility and cost, the possible ecological effects, and the social accept- 
ability of each countermeasure should be taken into account before it 
can be applied 7 ' 22 " 25 . In reality, the availability of countermeasure 
options for forest areas is limited to the restriction of access and 
removal of contaminated materials, and both of these options have 
their advantages and disadvantages 7,22 " 25 . Therefore, a combination 
of the above two options (zoning) could provide a realistic and 
balanced approach to the management of contaminated forests. 
Data revealing the radiocesium dynamics in Japanese forests are 
limited. At the time of writing, less than three years after the 
Fukushima accident, data are not necessarily sufficient fully to para- 
meterise models of radionuclide dynamic for longer-term predic- 
tions. However, the model used in this study was widely tested 
against data observed in forests contaminated by the Chernobyl 
accident, and we carefully parameterised and tested the model using 
data observed in the Fukushima area as fully as possible at this stage. 
Moreover, to develop effective countermeasures at an early stage, 
anticipation of the future dynamics of radiocesium in forest ecosys- 
tems is urgently required. We acknowledge that continued monitor- 
ing is necessary, together with iterative model parameterisation and 
validation in future years. Nevertheless, our simulation results at this 
early stage will assist in understanding and anticipating the future 
dynamics of radiocesium in Japanese forests. 

Methods 

The RIFE1 {Radionuclides in Forest Ecosystems) model was used in this study. The 
RIFE1 model was developed by Shaw and Belli 8 ' 1618 after the Chernobyl accident to 
track the radionuclide dynamics in forest ecosystems and has been tested and cali- 
brated using the data from eight forest sites in Europe 18 ' 9 . RIFE1 is one of the models 
tested in a model inter- comparison programme of International Atomic Energy 
Agency 1619 . RIFE1 is a relatively simple dynamic compartment model of radionuclide 
redistribution following a single deposition event and is designed to provide forecasts 
of the radionuclide dynamics in forest ecosystems to help develop decision making 
strategies following forest contamination (Supplementary Fig. Si). The model con- 
sists of five compartments {tree external, tree internal, litter, organic soil, and mineral 
soil), and these pools are linked by radionuclide fluxes {Supplementary Fig. Si). Each 
flux is described as a first-order rate process, which can be represented by a time- 
invariant rate coefficient. The model equations are as follows: 

dQi/df=Ii-(*i + fci+A)Qi 
dQ 2 / df = fei Qi + h Q4 - (*a + *) Q2 
dQ 3 /dt = I 2 + k 4 Q, + k 2 Q 2 - (k 5 + X) Q 3 
dQ 4 /df - k 5 Q$ - {h+h + X)Q 4 

dQ 5 /dt=hQ4-(k 7 +X)Q 5 

where Q l5 Q 2 , Q 3 , Q 4 , and Q 5 are the activity inventories of a radionuclide within the 
tree external, tree internal, litter (fresh dead leaves on the soil surface), soil surface 
organic soil (decomposed organic matter on the soil surface; F and H layers com- 
bined), and mineral soil components, respectively (Bq m" 2 ). fcj is the rate coefficient of 
radionuclide transfer between compartments (yr" ! ), and X is the rate constant of 
radioactive decay (yr -1 ) for the radionuclide being considered, t is the time (yr), and Ii 



and I 2 are the inputs of radionuclides to the tree external and litter components, 
respectively. The residence half-lives (Th, yr) used in this study are shown in 
Supplementary Table SI. The 137 Cs and l34 Cs dynamics were simulated in this study. 
The half-lives were converted into the appropriate rate coefficients using the fol- 
lowing relationship. 

k(orX)=ln(2)/T h 

Because it is difficult to estimate the leaf absorption parameter (hi), the rate was set 
to zero, as in the original study. The values for k 5 , k 6 , and k 7 were multiplied by the 
thickness of the layer when used. 

We used the inventories of radiocesium ( 137 Cs, 134 Cs) measured in September 201 1 
in Otama village (37 U 34'N, 140°18'E) in the Fukushima prefecture as the initial 
conditions for this simulation 26 (the location of the FDNPP is shown in 
Supplementary Fig. S2). The two most common forest types were included: evergreen 
needleleaf forest (mainly Cryptomeria japonica) and deciduous broadleaf forest 
(mainly Quercus serrata). The observed l37 Cs and l34 Cs distributions as percentages of 
the total inventory (i.e., 137 Cs + 134 Cs) were used as the initial conditions 
(Supplementary Table S2); hence, we simulated the radiocesium redistribution on a 
percentage basis. The major difference in the radiocesium inventory between ever- 
green needleleaf forests and deciduous broadleaf forests was that the contribution of 
the litter component was low for evergreen needleleaf forests and high for deciduous 
broadleaf forest, most likely because the trees in the deciduous forests were at a stage 
prior to leaf growth when the radioactive materials were discharged. The depth of the 
litter and organic soil for k 5 and k 6 were set to 0.017 and 0.012 m (N = 88), 
respectively, which were measured in Fukushima in the national soil carbon 
inventory project. In the field measurement data obtained, the portions of radio- 
cesium in the litter and organic soil were not distinguished; we assumed that 20% of 
the observed radiocesium was retained in the litter (L layer) and that 80% was in the 
organic soil (F and H layers). Because our field observation monitors the mineral soil 
layer over a depth of 0-0.2 m, we set the mineral soil layer thickness at 0.2 m for k 7 . In 
this study, we simulated the dynamics of 137 Cs, 134 Cs, and the total radiocesium in the 
tree components (Q L + Q 2 , Tree external and Tree internal in RIFE 1), soil surface 
organic layer (Q 3 + Q 4 , Litter and organic soil in the model), and mineral soil (Q 5 , 
Mineral soil in the model). The time courses of each raw component from the model 
are shown in Supplementary Figure S3. The transfers to understory vegetation and 
mushrooms were calculated using transfer factors in the RIFE1 model, but these 
transfers were not included in this simulation. Monthly time steps were adopted in the 
simulation. 

The model was parameterised using the data observed in the first and second years 
after the Fukushima accident 26 . A Bayesian calibration scheme, which is an optim- 
ization scheme that uses Monte Carlo sampling, was used to estimate the best para- 
meter sets (k 2 -k 7 ) and their uncertainty ranges (Supplementary Table Si). Details 
about the Bayesian calibration scheme are described elsewhere 27 . We conducted 
100000 iterations of sampling, and discarded the first 20000 iterations as burn-in to 
avoid the influence of the initial values. Using the original parameter set (geometric 
means of eight sites in Europe) the model overestimated the % of total radiocesium 
inventory in the soil surface organic layer and underestimated that in the mineral soil 
(Fig. S4). After the parameterisation, however, the model accurately reproduced the 
distribution of radiocesium in forest ecosystems one year after the fallout (Fig. S4). 
Because of the limited data, the parameters have large uncertainty ranges. 

The simulated results were validated by comparing estimated percentages of total 
radiocesium inventory in the second year with observed data at four needleleaf forest 
sites (Fig. S5 and Table S3) (ref. 26). The data used for validation were not used for 
model parameterisation. The inventories of radiocesium in tree external (E li2)3)4 ), tree 
internal (11,2,3,4), soil surface organic layer (Oi )2)3i4 ), and mineral soil (S 12)3i4 ) com- 
ponents were compared. The data were obtained in Kawauchi village (Ej, Ii, Oi, Si, 
and E 4 , 1 4 , 0 4 , S 4 ), Otama village (E 2 , 1 2 , 0 2 , S 2 ), and Tadami town (E 3 , 1 3 , 0 3 , S 3 ) 
(Table S3) (ref. 26). Overall, the simulation results were well correlated with observed 
data, showing high radiocesium inventories in mineral soil components and low 
inventories in tree and soil surface organic layers. However, relatively larger dis- 
crepancies between simulated and observed values were found for site 4 (E 4> 0 4 , S 4 ), 
where larger leaf and branch biomass and a thicker soil surface organic layer than 
other sites were reported 26 . These validations suggest that overall dynamics of 
radiocesium in Japanese forests were similar, and the RIFE1 model can accurately 
capture the overall dynamics of radiocesium in Japanese forests. However, since the 
dynamics of radiocesium in forests differ depending on forest conditions (e.g. tree 
biomass and soil surface organic layer), specific parameterisation would be needed for 
forests with distinctly different characteristics from the 'typical' forests of the 
Fukushima region. Also we emphasise that monitoring, together with iterative model 
parameterisation and validation, must be continued in future years. 

We combined the simulation results with a spatial radiocesium distribution map 
that was measured by an air-borne radiometric survey 1,2 . The map we obtained did 
not include the area near the FDNPP. We converted the map into a grid of 1 km 
squares and then extracted data from the areas containing forest. The evergreen 
needleleaf forest and deciduous broadleaf forests were distinguished using a 
vegetation database 28 . We multiplied the total radiocesium in a grid, as estimated by 
the air-borne survey (T RC , kBq m" 2 ), with our simulation results (e.g., I34 Q 5) ( _ 3 , the 
134 Cs in the mineral soil component 3 years after the fallout, as a proportion of the 
total initial deposition,%) and estimated the radiocesium dynamics in each grid. For 
example, the 134 Cs inventory in the mineral soil compartment 3 years after the fallout 
(kBq in" 2 ) can be calculated as T RC X 134 Q 5 , t = 3 /100 (e.g., 10.0 kBqm" 2 X 30%/100 
-3.0 kBq m" 2 ). 
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The dominant forest types in the contaminated area are evergreen needleleaf 
forests and deciduous broadleaf forests. The flora in Japanese forests is rich and 
diverse. The dominant tree species in evergreen needleleaf forests are Cryptomeria 
japonica (Japanese cedar) and Chamaecyparis obtusa (Japanese cypress), and those in 
deciduous broadleaf forests are Fagus japonica, Fagus crenata, Quercus serrata, and 
Quercus crispula. Brown forest soils (Cambisols and Andosols in the classification of 
the Food and Agriculture Organization) are distributed most widely in the area 
(approximately 70%), and Black soils (Andosols) are also present (approximately 
20%). 
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